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An effective characterization of chaotic conservative Hamihonian systems in terms of the curvature 
associated with a Riemannian metric tensor derived from the structure of the Hamiltonian has been 
extended to a wide class of potential models of standard form through definition of a conformal 
metric. The geodesic equations reproduce the Hamilton equations of the original potential model 
through an inverse map in the tangent space. The second covariant derivative of the geodesic 
deviation in this space generates a dynamical curvature, resulting in (energy dependent) criteria for 
unstable behavior different from the usual Lyapunov criteria. We show here that this criterion can 
be constructively used to modify locally the potential of a chaotic Hamiltonian model in such a way 
that stable motion is achieved. Since our criterion for instability is local in coordinate space, these 
results provide a new and minimal method for achieving control of a chaotic system. 
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I. INTRODUCTION 

The issue of controlling chaos has attracted great in- 
terest in the last two decades; much work has been done 
in the case of dissipative systems. Ott, Grebogi and York 
(OGY) pj developed a method by which chaos can be sup- 
pressed by making small time-dependent perturbations 
in order to shadow one of the infinitely many periodic 
orbits embedded in the chaotic attractor. This method 
has been extended by many following studies 0-11] and 
successfully applied to experimental systems @. 

As a result of the absence of attractors in conservative 
systems and the complex nature of the phase space com- 
bining regular and irregular regions, constructing an ef- 
fective method has remained a challenging question. Lai 
Ding and Grebogi (LDG) extended the OGY method to 
Hamiltonian systems by incorporating the notion of sta- 
ble and unstable directions at each periodic point. 

In relation to the above question, Chandre et al. 0- 
1^ have developed a method for constructing barriers in 
phase space to subdue chaotic behavior for large times. 
This method is based on the introduction of a specifically 
designed small control term which changes the dynam- 
ics from chaotic to regular behavior. Zhang et al. Q 
have introduced a method for controlling chaos in two- 
dimensional Hamiltonian system which they called adap- 
tive integrable mode coupling based on the separation of 
the system into two coupled subsystems one of which is 
stable and the other unstable; when the unstable system 
comes into the vicinity of the integrable system the con- 
ditions are reset resulting in an effective adaptive con- 
trol. Cartwright et al. [IQ] have constructed a method 
which permits stabilization of KAM islands through for- 
ward iteration of the orbits, and transforming them into 
global attractors of the embedded system. Ciraolo et 
al. [HI have used a method of adding a small pertur- 
bation to modify systems, such as nonlinear plasmas, to 



follow more regular motion. Additional efficient methods 
have been proposed in [T^ - [T6j : some of these methods 
require tracking the trajectories while others involve the 
addition of interacting terms for which the criteria are 
somewhat sensitive. 

Geometric approaches for the analysis of the stability 
of a given Hamiltonian have been widely discussed [17|- 
[29| . for which curvature of a manifold is associated with 
stability. 

The pioneering work of Oloumi and Teychenne [s^l 
proposes a stabilization method by considering the Gaus- 
sian curvature of the potential energy as a source of 
chaos. Even if the condition of negative curvature of the 
potential and instability of the dynamics are not com- 
pletely equivalent they were able to successfully control 
the instability by avoidance of negative curvature regions 
of the potential energy (ANCRP). 

In this work we make use of a recently developed ge- 
ometrical criterion [25| which is highly sensitive to in- 
stability in Hamiltonian systems, and has been shown 
to be in complete agreement with the results of the nu- 
merical technique of surface of section (Poincare plot) 
and more effective than other geometric methods [27|. 
This criterion is based on an equivalence between motions 
generated by Hamiltonian in the standard form with 
quadratic kinematic terms and additive potentials (which 
we shall call the Hamilton description) and a Hamiltonian 
in which the dynamics is described by a metric- type func- 
tion of coordinates multiplying the momenta in bilinear 
form (which we shall call the Gutzwiller description) (sTI - 
Issj . The mapping between these equivalent descriptions 
was first introduced by Appel 36-39]. The criterion for 
stability developed in 25] establishes an inverse map- 
ping of the motion described completely geometrically in 
the Gutzwiller space back to the Hamilton description, 
carrying with it the covariance under diffeomorphisms 
that is a property of the Gutwiller dynamics. The result- 
ing orbits in the Hamilton description, with the special 
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choice of coordinates serving as the basis for the Appel 
type relation, reduce to the standard description of the 
Hamilton orbits under the Hamilton equations. In gen- 
eral geometric form, as geodesic equations, however, they 
are subject to analysis in terms of geodesic deviation, and 
the resulting formula can be reduced to a computation in 
the special coordinates of a new (symmetric) matrix val- 
ued criterion for stability [25[ which has been shown to 
be very effective in a wide range of examples [13, [13, 113 ■ 



In the analysis of these systems, it was shown that the 
presence of negative eigenvalues for the stability matrix 
in the admissible physical region (for which E > V) re- 
sults in a chaotic type Poincare plot. Since the condition 
is local in coordinate space, the result necessarily implies 
some degree of ergodicity. There has, however, been no 
proof that the existence of negative eigenvalues is not 
strongly associated with the structure of the dynamics 
more globally, as is often the case, for example, for an- 
alytic function theory, where the presence of a complex 
pole implies a distortion of a relatively large region of 
the complex plane. The results of our application to the 
control of dynamical systems, however, indicates that the 
regions of negative eigenvalue correspond to a very local- 
ized phenomena. The observed instability they induce 
appears to be due to the passage of orbits through these 
regions, and they are not strongly correlated with the 
more global structure of the dynamics. 



What we have done here is to identify the regions of 
negative eigenvalue in the coordinate space, and mod- 
ify the Hamiltonian in these regions locally, removing 
the source of instability either by removing the insta- 
bility inducing terms or varying the coupling to these 
terms to bring the Hamiltonian closer to an integrable 
form in these regions. The remaining part of the space 
is left to develop according to the full, nonlinear, sym- 
metry breaking, evolution. The effect on the Poincare 
plots, as we show here, is dramatic. This result provides 
strong evidence that chaotic Hamiltonian systems derive 
their properties from sources that may be thought of as 
highly localized in the regions of negative eigenvalues for 
the stability matrix. In addition to giving an interest- 
ing insight into the nature of Hamiltonian chaos, it also 
provides an effective control method in which the Hamil- 
tonian dynamics is left to exercise its full function in a 
large domain of configurations, but is only constrained in 
prescribed local regions of the configuration space. 



In the next section, we review the basic ideas underly- 
ing the geometrical criterion, and in the following section 
we give numerical results for the application of our con- 
trol procedure for the cases of an oscillator with broken 
symmetry and a potential which may be derived from the 
Toda form. 



II. THEORY 

It has been shown [25|, [s^ that a Hamiltonian system 
of the form (we use the summation convention) 



H = 



2AI 



Vix) 



(1) 



where y is a function of space variables alone, can be put 
into the equivalent form 



Hg = 



(2) 



where gij is conformal and is a function of the coordinates 
alone. One can easily see that the orbits described by the 
Hamilton equations for ([2|) coincide with the geodesies on 
a Riemannian space associated with the metric gij, i.e., it 
follows directly from the Hamilton equations associated 
with ^ that [Sli] 
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where the connection form F™" is given by 
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and is the inverse of gij . 
For a metric of conformal form 

gij = ^Sij , 



(3) 



(4) 



(5) 



with inverse g^^ = ^p~^5ij, O'n. the hypersurface defined by 
Hq = H = E = constant, and assuming that the observ- 
able momenta are the same at every t, the requirement 
of equivalence implies that [25l | 



E 



E~V{x)' 



(6) 



To see that the Hamilton equations obtained from ([2]) 
can, be put into correspondence with those obtained from 
the Hamiltonian of the potential model ([1]), we first note, 
from the Hamilton equations for ([2]), that 



• ^9Hg _ 1 
dp' M 

We then define the velocity field 
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coinciding formally with one of the Hamilton equations 
implied by ([T]) , for which we label the coordinates x^ . 

To complete our correspondence with the dynamics in- 
duced by ([T]) , consider the Hamilton equation for p' gen- 
erated by Hg, 
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FIG. 1: Effect of control on the model of Eq. (|19[) . The first column shows the physical region (closed black curves) corresponding 
to the sequence of energies, corresponding to a, b, c, d. The regions of negative eigenvalues (instability) are shown in 

gray. The second column shows the Poincare plots for the uncontrolled system. The third column shows the Poincare plots for 
the controlled system (Eq. (|2ip with r = 0.4 for all cases). The (red) dashed line in column one shows the boundary of the 
control modification. 



With 



where 



this becomes 



(10) 



(11) 



Eq. ([TU|) has the form of a geodesic equation, with 
a truncated connection form |25| . Note that perform- 
ing parallel transport on the local flat tangent space of 
the Gutzwiller manifold (for which r™" and gij are com- 
patible), the resulting connection, after raising the ten- 



sor index to reach the Hamilton manifold, is exactly the 
"truncated" connection ([TT|) . 

Substituting (P and ^ into ^ and (dJ), the Kro- 
necker deltas identify the indices of i™ and i"; the re- 
sulting square of the velocity cancels a factor of {E — 
V)~^ , leaving the Hamilton- Newton law derived from 
the Hamilton equations directly from ([1]). Eq. ((TU)) is 
therefore a geometrically covariant form of the Hamilton- 
Newton law, exhibiting what can be considered an under- 
lying geometry of standard Hamiltonian motion. 

Since the coefficients constitute a connection 

form, they can be used to construct a covariant deriva- 
tive. It is this covariant derivative which must be used to 
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FIG. 2: Effect of control on the model of Eq. (|20|) . The first column shows the physical region (closed black curves) corresponding 
to the sequence of energies, 0.5, 0.6, 0.8, 1.0, 1.5 corresponding to a, 6, c, d, e. The regions of negative eigenvalues (instability) are 
shown in gray. The second column shows the Poincare plots for the uncontrolled system. The third column shows the Poincare 
plots for the controlled system (Eq. (|22p with a = 0.5 for all cases). The (red) dashed line in column one shows the boundary 
of the control modification. 
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compute the rate of transport of the geodesic deviation 
_ ^it _ along the (approximately common) motion 

of neighboring orbits in the Hamilton manifold, since it 

follows the geometrical structure of the geodesic curves 

on {x^}. 

The relation 

e = -2Mi^i^C - (12) 



obtained from dTOt . can be factorized in terms of the 
covariant derivative 



One obtains 



D 



M 



DMt^ 
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(13) 



(14) 



where the index M refers to the connection (ITT|) . and 
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dM: 
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(15) 

corresponds to the curvature associated with the connec- 
tion form This expression does not coincide with 
the curvature of the Gutzwiller manifold (given by this 
formula with in place of Afg„j), but is a dynami- 
cal curvature which is appropriate for geodesic motion in 

With the conformal metric in noncovariant form (O,® 
(in the coordinate system in which ^ is defined) , the dy- 
namical curvature (1151) can be written in terms of deriva- 
tives of the potential V, and the geodesic deviation equa- 
tion becomes 



where the matrix V is given by 
3 dV dV 



and 



1 d'^V 
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(16) 



(17) 



(18) 



with = x^ , defining a projection into a direction or- 
thogonal to . 

Instability should occur if at least one of the eigenval- 
ues of PVP is negative, in terms of the second covariant 
derivatives of the transverse component of the geodesic 
deviation. This condition is easily seen to be equivalent 
to the same condition imposed on the spectrum of V, and 
is thus independent of the direction of the motion on the 
orbit. 

The condition implied by the geodesic deviation equa- 
tion (jl6p . in terms of covariant derivatives, in which the 



orbits are viewed geometrically as geodesic motion, is a 
new condition for instability |25| . based on the underly- 
ing geometry, for a Hamiltonian system of the form ([T|), 
providing new insight into the structure of the unstable 
and chaotic behavior of Hamiltonian dynamical systems. 
The method has been successfully applied to many po- 
tential models [13, US 113, EO] ■ We now apply the method, 
which is simple and straightforward, since the instability 
criterion is defined locally in coordinate space, to extract 
from a chaotic Hamiltonian the part of the dynamics in 
the regions for which V has negative eigenvalues which 
lead to unstable motion . We show by direct simulation 
that this procedure is remarkably effective. 

The results confirm the notion of the locality of our cri- 
terion, since the (local) removal of parts of the Hamilto- 
nian inducing instability by this criterion strongly affects 
the global stability of the motion. 



III. RESULTS AND DISCUSSION 

In the following we give some examples of Hamiltonian 
motion which are unstable and exhibit chaotic behavior. 
In two dimensions, the terms in the potential that break 
rotational symmetry have a form, in these examples, that 
induces chaotic behavior. The matrix V computed over 
the physically accessible region exhibits certain regions 
with negative eigenvalues. We have arranged our simu- 
lation to replace the Hamiltonian in the regions of nega- 
tive eigenvalues by a Hamiltonian in which the symmetry 
breaking terms are either removed or for which the cou- 
pling coefficient is decreased, but in regions of positive 
eigenvalues, the Hamiltonian remains in its original form. 
The changes are carried in relatively small regions of the 
configuration space in the neighborhood of the regions of 
negative eigenvalues. 

We take for illustration here a simple and important 
case of coupled harmonic oscillators with perturbation 



V{x,y) = i(a;2 + y2) + 6xV 
and a generalization of the Toda potential 



(19) 



v{x, y) + y') + ^'y - ly' + + \y'. (20) 

The first of these is known to generate chaotic behav- 
ior; as we have shown previously, the transition to chaotic 
behavior as a function of energy and of the coupling to 
the perturbation is well described by our geometric cri- 
terion 25]. In the present work, we set the coupling to 
the chaos inducing perturbation to zero in the regions for 
which the geometric criterion results in negative eigenval- 
ues. 

In our previous study of the second example Eq. ((20)) 
we studied the sensitivity of the dynamical behavior to 
the choice of energy [23| . In the present work, at the en- 
ergies for which the system exhibits chaotic behavior we 
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have removed the chaos inducing perturbation in the re- 
gions for which the geometric criterion resuhs in negative 
eigenvalues. 

Fig. 1 shows the effects on the dynamical behavior of 
the change of the coupling to a stable value in the re- 
gions of negative eigenvalues for the perturbed oscillator 
potential (IT^ . The results for different values of energy 
are shown; these values correspond to different sizes of 
the physically accessible region. Note that the radius of 
the region of positive eigenvalues does not change in this 
example. 

The first row shows the physical region and the lo- 
cation of negative eigenvalues. The interior of the cir- 
cle corresponds to a region in which no negative eigen- 
value occur, chosen in this way to simplify the compu- 
tation. The second row shows the surface of section 
Poincare plot (surface of section) of the original uncon- 
trolled Hamiltonian indicating chaotic dynamics. The 
third row shows the Poincare plot generated by the con- 
trolled potential: 



V{x,y) = 




where r stands for the radius of the region of positive 
eigenvalues. 

Fig. 2 shows the effect of control on second system 
(1^(71) , using the controlled potential: 

(22) 



where a stands for the limit (a horizontal dashed line) of 
the region of positive eigenvalues. 

In both cases, for sufficiently high energies the uncon- 
trolled systems become less stable and a chaotic signature 
appears. We examine the method for different values of 
the energy, all corresponding to chaotic motion of the un- 
controlled system. One can easily see that the Poincare 
plots in controlled systems present almost completely reg- 
ular motion. 



IV. CONCLUSION 



We see from these computations that in regions of pos- 
itive eigenvalues the Hamiltonian with chaos inducing 
terms do not generate instability. If the chaos induc- 
ing terms are removed or decreased in coupling in the 
regions of negative eigenvalues, the chaotic motion of the 
system is reduced or disappears. We conclude that the 
chaotic behavior of the system is associated with local 
properties of the Hamiltonian, and that our local crite- 
rion is effective in identifying these regions. Furthermore, 
our procedure provides an effective method for control of 
chaotic systems of this type. 
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